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Abstract 

This work discusses the perspectives to observe fluxes of high energy astrophysical neutrinos with the planned 
km 3 telescopes. On the basis of the observations of GeV and TeV 7-rays, and of ultra high energy cosmic rays, 
it is possible to construct well motivated predictions that indicate that the discovery of such fluxes is probable. 
However the range of these predictions is broad, and the very important opening of the "neutrino window" on 
the high energy universe is not guaranteed with the current design of the detectors. The problem of enlarging the 
detector acceptance using the same (water/ice Cherenkov) or alternative (acoustic/radio) techniques is therefore 
of central importance. 



1. Introduction 

The fundamental significance of the opening of 
the new "window" of high energy v astronomy 1 
for the observation of the universe is beyond dis- 
cussion. Neutrinos have properties that are pro- 
foundly different from those of photons, and ob- 
servations with this new "messenger" will allow us 
to develop a deeper understanding of known as- 
trophysical objects, and also likely lead to the dis- 
covery of new unexpected classes of sources, in the 
same way as the opening of each new window in 
photon wavelength has lead to remarkably inter- 
esting discoveries. 

If there are few doubts that in the future neu- 
trino astronomy will mature into an essential field 
of observational astrophysics, it is less clear how 
long and difficult the history of this development 



1 Neutrino Astronomy for solar and gravitational collapse 
SuperNovae v has already been remarkably successful. This 
discussion will only consider high energy (E L 
neutrinos. 
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will be, and in particular what will be the scientific 
significance of the results of the planned km 3 tele- 
scopes. These telescopes are "discovery" instru- 
ments and only a posteriori, when their data analy- 
sis is completed, we will be able to appreciate their 
scientific importance. It is however interesting to 
attempt an estimate, based on our present knowl- 
edge, of the expected event rates and number of 
detectable sources. 

As a warning, it can be amusing to recall that 
when (in june 1962) Bruno Rossi, Riccardo Gi- 
acconi and their colleagues flew the first X-ray 
telescope [1,2] opening the X-ray photon window 
to observations, the most promising X-ray source 
was the sun, followed by the moon (that could scat- 
ter the solar wind) . It is now known that the moon 
surface does indeed emit X-rays, but in fact it ap- 
pears as a dark shadow in the X-ray sky, because 
it eclipses the emission of the ensemble of Active 
Galactic Nuclei (AGN) that can now be detected 
with a density of ~ 7000 per square degree [3] . So 
the X-ray sky was very generous to the observers, 
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and embarassed the theorists who had worked on 
a priori predictions. 

The lesson here is that predictions about the un- 
known are difficult. One should expect the unex- 
pected, and it is often wise to take (calculated) sci- 
entific risks. We all hope that history will repeat 
itself, and that also the v sky will be generous to 
the brave scientists who, with great effort, are con- 
structing the new telescopes. 

On the other hand, we are in the position to make 
some rather well motivated predictions about the 
intensity of the astrophysical v sources, because 
high energy v production is intimately related to 
7-ray emission and cosmic ray (c.r.) production, 
and the observations of cosmic rays and high en- 
ergy (GeV and TeV range) photons do give us very 
important guidance for the prediction of the neu- 
trino fluxes. 

This review will concentrate only on high energy 
astrophysics. There are other important scientific 
topics about km 3 telescopes that will not be cov- 
ered here. A subject of comparable importance is 
"indirect" search for cosmological Dark Matter via 
the observations of high energy v produced at the 
center of the sun and the earth by the annihila- 
tion of DM particles [4] . The km 3 telescopes can 
also be used to look for different "New Physics" 
effects such as the existence of Large Extra dimen- 
sions [5] . The instruments also have the potential 
to perform interesting interdisciplinary studies. 

2. Components of the Neutrino Flux 

The observable neutrino flux can be schemati- 
cally written as the sum of several components: 
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The first two components describe atmospheric 
neutrinos that are generated in cosmic ray showers 



in the Earth's atmosphere. They are an important 
foreground to the more interesting observations 
of the astrophysical components. Atmospheric i/'s 
can be split into two components, the first "stan- 
dard" one is due to the decay of charged pions and 
kaons, while the second one is due to the weak 
decays of short lived (hence "prompt") particles 
containing heavy quarks, with charmed parti- 
cles accounting for essentially all of the flux. The 
prompt contribution is expected to be dominant 
in the atmospheric neutrino fluxes at high energy. 
This component has not yet been identified, and 
its prediction is significantly more difficult than 
the standard flux, because of our poor knowledge 
of the dynamics of charmed particles production 
in hadronic interactions. Figure 1 shows the angle 
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Fig. 1. Neutrino fluxes. The solid lines labeled u e , f M and 
v T show the angle averaged fluxes of atmospheric neutrinos 
and anti-neutrinos as a function of energy. The associated 
dashed lines are the fluxes neglecting the "prompt" con- 
tribution of charmed hadrons. The lines labeled i^-GZK 
describes GZK neutrinos [6]. The points (and fitted lines) 
describe CR direct and indirect measurements. The hori- 
zontal lines labeled Amanda and Icecube are the current 
limit and the predicted sensitivity of km 3 v telescopes for 
an isotropic v flux. The line labeled Anita-Lite is the limit 
on the v flux obtained with radio methods [7] . The isotropic 
7— ray flux measured by Egret [8] is also shown. 

averaged atmospheric v fluxes [9] with an estimate 
of the "prompt" component [10] that "overtakes" 
the standard one at E v ~ 10 TeV for v e and 
- 100 TeV for v^. 

The next two components of the neutrino fluxes 
are "diffuse" fluxes coming from v production in 
interstellar space in our own Galaxy, and in in- 
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tergalactic space. The diffuse Galactic emission is 
due to the interaction of cosmic rays (confined in- 
side the Milky Way by the galactic magnetic fields) 
with the gas present in interstellar space. The an- 
gular distribution of this emission is expected to 
be concentrated in the galactic plane, very likely 
with a distribution similar to the one observed for 
GeV photons by EGRET [11]. The extragalactic 
diffuse emission is dominated by the decay of pi- 
ons created in p-y interactions by Ultra High En- 
ergy protons (E p > 6 x 10 19 eV) interacting with 
the 2.7K cosmic radiation. These "GZK" neutri- 
nos are present only at very high energy (with the 
flux peaking at E v ~ fO 18 eV). 

Together with the diffuse fluxes one expects the 
contribution of an ensemble of point-like (or quasi- 
point like) sources of galactic and extragalactic ori- 
gin. Neutrinos travel along straight lines and al- 
low the imaging of these sources. It is expected 
that most of the extragalactic sources will not be 
resolved, and therefore the ensemble of the extra- 
galactic point sources (with the exception of the 
closest and brightest sources) will appear as a dif- 
fuse, isotropic flux that can in principle be sepa- 
rated from the atmospheric v foreground because 
of a different energy spectrum, and flavor compo- 
sition. 



3. Neutrinos, Photons and Cosmic Rays 

In the standard mechanism for the production 
of high energy v in astrophysical sources a popu- 
lation of relativistic hadrons (that is cosmic rays) 
interacts with a target (gas or radiation fields) cre- 
ating weakly decaying particles (mostly n ± , and 
Kaons) that produce v in their decay. The energy 
spectrum of the produced neutrinos obviously re- 
flects the spectrum of the parent cosmic rays. A 
well known consequence of the approximate Feyn- 
man scaling of the hadronic interaction inclusive 
cross sections is the fact that if the parent c.r. have 
a power law spectrum of form (f> CI ~ K cr E~ a , 
and their interaction probability is energy indepen- 
dent 2 the v spectrum, to a good approximation, 

2 In the most general case c.r. of different rigidity p/Z 
diffuse in different ways inside the source and have different 



is also a power law of the same slope. 

The current favored models for 1st order Fermi 
acceleration of charged hadrons near astrophysical 
shocks predict a generated spectrum with a slope 
a ~ 2. This expectation is in fact confirmed by the 
observations of young SNR by HESS, and leads the 
expectations that astrophysical v sources are also 
likely to have power law spectra with slope close to 
2. Such a slope is also predicted in Gamma Rays 
Bursts models [20] where relativistic hadrons inter- 
act with a power law photon field. It is important 
to note that the high energy cutoff of the parent 
c.r. distribution is reflected in a much more gradual 
steepening of the v spectrum for E v > 0.01 E Ptmax . 

The hadronic interactions that are the sources 
of astrophysical neutrinos also create a large num- 
ber of 7r° and rj particles that decay in a 77 mode 
generating high energy 7-rays. In general, it is pos- 
sible that these photons are absorbed inside the 
source, and the energy associate with them can 
emerge at lower frequency, however for a transpar- 
ent source the relation between the photon and 
neutrino fluxes is remarkably robust. For a power 
law c.r. spectrum, the 7-rays are also created with 
a power law spectrum of the same slope. The ap- 
proximately constant Z//7 ratio is shown in fig. 2 
as a function of the slope a. The vjy ratio is ap- 
proximately constant in energy when E is much 
smaller than the c.r. energy cutoff, and depends 
only weakly on the hadronic interaction model. 

At the source, the relative importance of the dif- 
ferent v types is: 

{v e ,v e , v^V^ v T ,V T } ~ {l + e,l-e,2,2,0,0} 

where e ~ 0.1 depends on the slope of the spec- 
trum, the relative importance of protons and heavy 
nuclei in the c.r. parent population, and the nature 
of the c.r. target (gas or radiation field). The ra- 
tio v^/fe — 2 at the source is a well known conse- 
quence of the fact that the chain decay of a charged 
pion: 7r + — > + followed by fi + — > e + + ^ e +F M 
generates two v of ^(-flavor and one of e-type. The 
presence of a v ^ and an in the final state insures 
that the ratio v^jv^ ~ 1, while the ratio v e jv e is 
controled by the relative importance of 7r + over 7r~ 



space distributions. For a non homogeneous target, this can 
be reflected in an energy dependent interaction probability. 
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Fig. 2. Ratio <j> v /<j>y between the neutrino and photon fluxes 
at E nj 3 TeV, for a power law c.r. population of slope 
a and p/n ratio of 10% interacting with a low density 
gas target. Photon absorption is neglected. The calculation 
uses the Sibyll [12] hadronic interaction model. The thick 
line sums over all neutrino type. The v flavor is calculated 
at the source before the inclusion of oscillations. 

production that is not symmetric for a p rich c.r. 
population. 



with respect to the typical size of astrophysical 
sources. Therefore oscillations are negligible for at- 
mospheric v above E v ~ 1 TeV, but can be safely 
averaged for (essentially all) astrophysical neutri- 
nos. After space averaging, the oscillation proba- 
bility matrix can be written as: 



(P( Va - V )) = (P(V a - V )) = J2 \Uaj\ 2 \U t 



0j\ 




(1) 



where we have used the best fit choice for the mix- 
ing parameters (6> 23 = 45°, 9 12 ~ 32.3°, 13 = 0). 
The most important consequence of (1) is a ro- 
bust prediction, valid for essentially all astrophys- 
ical sources, for the flavor composition of the ob- 
servable v signal: 

{v e + V e , Uft + V f _ l , v T + V T } = {1, 1, 1} 



5. The Gamma ray sky 



4. Neutrino Oscillations 

Measurements of solar and atmospheric neu- 
trinos have recently established the existence of 
neutrino oscillations, a quantum-mechanical phe- 
nomenon that is a consequence of the non-identity 
of the v flavor {v ei v^, v T } and mass {v\, vz\ 
eigenstates. A v created with energy E and fla- 
vor a can be detected after a distance L with a 
different flavor (3 with a probability that depends 
periodically on the ratio L/E. This probabil- 
ity oscillates according to three frequencies that 
are proportional to the difference between the v 
squared masses Am 2 fe , and different amplitudes 
related to the mixing matrix U a j that relates the 
flavor and mass eigenstates. The shortest (longest) 
oscillation length, corresponding to the largest 
(smallest) |Am 2 |) can be written as: 

A12 = (4nE u )/\Aml 2 \ ~3.1x 10 12 £ T eV cm, 
A 23 = (4nE v )/\Aml 3 \ ~ 0.99 x 10 11 .BtcV cm. 

These lengths are long with respect to the Earth's 
radius (Rq ~ 6.371 x 10 8 cm), but are very short 



The TeV energy range (the highest avaliable) for 
photons is the most interesting one for neutrino as- 
tronomy, because the v signal for the future tele- 
scopes is expected to be dominated by neutrinos of 
one (or a few) order(s) of magnitude higher energy. 
Recently the new Cherenkov 7-ray telescopes have 
obtained remarkable results, and the catalogue of 
high energy gamma ray sources has dramatically 
increased. Of particular importance has been the 
scan of the galactic plane performed by the HESS 
telescope [13,14], because for the first time a cru- 
cially important region of the sky has been ob- 
served with an approximately uniform sensitivity 
with TeV photons. 

The three brightest galactic TeV sources de- 
tected by the HESS telescope have integrated 
fluxes above 1 TeV (in units of 10~ n (cm 2 s) _1 ) 
of approximately 2.1 (CRAB Nebula), 2.0 (RX 
J1713.7-3946) and 1.9 (Vela Junior) 

The fundamental problem in the interpretation 
of the 7-ray sources is the fact that it is not known 
if the observed photons have hadronic (ir° decay) 
or leptonic (inverse Compton scattering of rela- 
tivistic electrons on radiation fields) origin. If the 
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leptonic mechanism is acting, the hadronic com- 
ponent is poorly (or not at all) constrained, and 
the v emission can be much smaller than the 7-ray 
flux. For the hadronic mechanism, the v flux is at 
least as large as the photon flux, and higher if the 
7-rays in the source are absorbed (see sec. 3). 

The 7-ray TeV sources, belong to several dif- 
ferent classes. The Crab is a Pulsar Wind Neb- 
ula, powered by the spin down by the central neu- 
tron star. The emission from these objects is com- 
monly attributed to leptonic processes, and in par- 
ticular the Crab is well described by the Self Syn- 
chrotron Comptons model (SSC). The next two 
brightest sources [15,16,17] are young SuperNova 
Remnants (SNR), and there are good reasons to 
believe that the photon emission from these ob- 
jects is of hadronic origin. The 7 spectra are power 
law with a slope a ~ 2.0-2.2 which is consistent 
with the expectation of the spectra of hadrons ac- 
celerated with 1st order Fermi mechanism by the 
SN blast wave. The extrapolation from the photon 
to the neutrino flux is then robust, the main uncer- 
tainty being the possible presence of a high energy 
cutoff in the spectrum. 

Other TeV sources that are very promising for 
v astronomy are the Galactic Center [18] with a 
measured flux ~ 0.2 (same units: 10 -11 (cm 2 s) -1 ) 
and the micro-Quasar LS5039 [19], the weakest 
observed TeV source with a flux <~ 0.12. These 
sources are not particularly bright in TeV 7-rays, 
but there arc reasons to believe that they could 
have significant internal absorption for photons, 
and therefore have strong v emission. 

In general, under the hypothesis of (i) hadronic 
emission and (ii) negligible 7 absorption (that cor- 
respond to <f> v ~ 7 ), even we assume that the 7 
and v spectra extends as a power law up to 100 TeV 
or more, the HESS sources are just at (or below) 
the level of sensitivity of the new v telescopes as 
will be discussed in the next section. 

For extra-galactic sources the constraints from 
TeV photon observations are less stringent because 
very high energy photons are severely absorbed 
over extra-galactic distances due to 77 — > e+e~ 
interactions on the infrared photons (from red- 
shifted starlight) that fill intcrgalactic space. Ac- 
tive Galactic Nuclei (AGN) are strongly variables 
emitters of high energy radiation. The EGRET de- 



tector on the CGRO satellite has detected over 90 
AGN of the Blazar class 3 with £ 7 > 100 MeV 
[31] . The brightest AGN sources in the present TeV 
7-ray catalogue, are the nearby (z ~ 0.03) AGN's 
Mkn 421 and 501 that are strongly variable on 
all the time scales considered from few minutes to 
years. Their average flux can be, during some pe- 
riods of time, several times the Crab. The extrap- 
olation to the v flux has significant uncertainties 
because the origin (hadronic or leptonic) of the de- 
tected 7 rays is not established. 

Leading candidates as sources of high energy 
neutrinos are also Gamma Ray Bursts (GRB) 
[20,21,22,23]. It is possible that individual GRB 
emits neutrino fluences that are sufficiently high 
to give detectable rates in the neutrino telescopes. 



6. Point Source Sensitivity 

The most promising technique for the detection 
of neutrino point sources is the detection of v— 
induced muons [p, j). These particles are produced 
in the charged current interactions of and 
in the matter below the detector. To illustrate this 
important point we can consider a "reference" v 
point source 4 with an unbroken power law spec- 
trum of slope a — 2.2, and an absolute normaliza- 
tion (summing over all v types) $„(> 1 TeV) = 
10~ n (cm 2 s) _1 . The reference source flux corre- 
sponds to approximately one half the flux of the 
two brightest SNR detected by HESS. The event 
rates from the reference source of v interactions 
with vertex in the detector volume, and for v— 
induced muons are shown in fig. 3, The event rates 



3 Blazars are AGN that emit one jet nearly parallel (within 
an angle T -1 , where T is the bulk Lorcntz factor of the 
jet) to the line of sight. 

4 In this work we have chosen to characterize the normal- 
ization of a v point source as the flux (summed over all v 
types) above a minimum energy of E v m ; n = 1 TeV. The 
reason for this choice is that it allows an immediate com- 
parison with the the sources measured by TeV 7— ray tele- 
scopes, that are commonly stated as flux above i? 7j m i n = 
1 TeV. In case of negligible 7 absorption one has <j>j ~ <j> v . 
Since we consider power law fluxes, it is trivial to restate 
the normalization in other forms. As discussed later the 
km 3 telescopes sensitivity peaks at E v ~ 20 TeV. 
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Fig. 3. Event rates from a v point source. The source has 
spectrum 4> v = K E^ a with slope 2.2 and normalization 
4>v(> 1 TeV) = 10 — 11 (cm 2 s)~ 1 equaly divided among 6 
v types. The thin lines describe the rate of iz-interactions 
in the detector volume for e, fi and r like CC interactions. 
The thick lines are fluxes of ^-induced muons with 1 GeV 
and f TeV threshold. 
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Fig. 4. Event rates for e— like CC interactions (thin lines) 
and v— induced muons with a 1 GeV threshold (thick lines) . 
The source is the same as in fig. 3). The solid (dashed) 
curve are for no (maximum) absorption in the Earth. 

for e, [i and r like v interactions are: 10.3, 9.6 
and 2.9 events/(km 3 yr) (assuming a water filled 
volume), while the ^-induced muon flux is cf>^ ~ 
5.6 (km 2 yr) -1 . The signal depends on the zenith 
angle of the source, because of v absorption in the 
Earth as shown in fig. 4. The rate of e-like events 
has a small contribution from the "Glashow reso- 
nance" (the process v e + er — > W~), visible as a 



peak at E v 



r/(2m e ) - 6 x 10 6 GeV. Neglect- 



ing the Earth absorption the resonance contributes 
a small rate ~ 0.07 (km 3 yr) -1 to the source signal, 



however this contribution quickly disappears when 
the source drops much below the horizon, because 
of absorption in the Earth. 

The n | rate is much easier to measure and to 
disentangle from the atmospheric v foreground. A 
crucial advantage is that the detected muon allows 
a high precision reconstruction of the v direction, 
because the angle 9^ v is small. 

In general the relation between the neutrino and 
the ^-induced muon fluxes is of order: 



(1-5) 
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The exact relation (shown 5 in fig. 5) depends on 
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Fig. 5. Flux of v— induced muons, calculated for a v point 
source with a power law spectrum of slope o, plotted as 
a function of the high energy cutoff of the parent proton 
spectrum. The thick (thin) lines are calculated for a thresh- 
old E u 

min — 1 GeV (1 TeV). The absolute normalization 
is chosen so that in the absence of cutoff the total v flux 
above f TeV has a flux $„ = IO" 11 (cm 2 s)" 1 . 

(i) the slope of the v spectrum, (ii) the presence 
of a high energy cutoff, (iii) the threshold energy 
used for fj, detection and (iv) (more weakly) on the 
zenith angle of the source (because of absorption 
effects) . 

The peak of the "response" curve for the v— 
induced muons (see fig. 3) is at E v ~ 10 TeV 
(20 TeV for a threshold of 1 TeV for the muons), 



5 Our numercal results are in excellent quantitative agree- 
ment with the calculation of Costantini and Vissani [24] 
for the case a = 2.2, when we consider a spectrum with 
no— cutoff. The calculation in [24] underestimates the effect 
of a cutoff in the proton parent spectrum. 
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with a total width extending approximately two 
orders of magnitude. In other words the planned 
v telescopes should be understood mostly as tele- 
scopes for ~ 10-100 TeV neutrino sources. This 
energy range is reasonably well "connected" to the 
observations of the atmospheric Cherenkov 7-ray 
observations that cover the 0.1-10 TeV range in 
-Ely, and the v sources observable in the planned 
km 3 telescopes are likely to be appear as bright 
objects for TeV 7-ray instruments. 

Because of the background of atmospheric 
muons, the ^-induced jit's are only detectable when 
the source is below the horizon. This reduces the 
sky coverage of a telescope as illustrated in fig. 6, 
that gives the time averaged signal obtained by 
two detectors, placed at the south-pole and in the 
Mediterranean sea, when the reference source we 
are discussing is placed at different celestial decli- 
nations. For a south pole detector the source re- 
mains at a fixed zenith angle cos0 ze mth = — sin (5, 
and the declination dependence of the rate is only 
caused by difference in v absorption in the Earth 
for different zenith angles. The other curves in- 
cludes the effect of the raising and setting of the 
source below the horizon. 

Some of the most promising galactic sources are 
only visible for a neutrino telescope in the northern 
Hemisphere since the Galactic Center is at declina- 
tion 5 ~ -29°. The interesting source RX J1713.7- 
3946 is at 6 = -39°, Vela Junior at S ~ -46°. 

6.1. Background Estimates 

Since the prediction of the signal size for from 
the (expected) brightest v sources is of only few 
events per year, it is clearly essential to reduce all 
sources of background to a very small level. This 
is a possible but remarkably difficult task. 

The background problem is illustrated in fig. 7 
that shows the energy spectrum of the muon sig- 
nal from our "reference" point source, comparing 
it with the spectrum of the atmospheric v fore- 
ground integrated in a small cone of semi-angle 
0.3°. The crucial point is that the energy spectrum 
from astrophysical sources is harder than the at- 
mospheric v one, with a median energy of approx- 
imately 1 TeV, an order of magnitude higher. It 
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Fig. 6. Average flux of up-going muons from a v point 
source of fixed luminosity, plotted as a function of decli- 
nation for two v telescopes placed at the south pole (thick 
line) and in the mediterranean sea (thin line). The flux is 
considered detectable only when the source is below the 
horizon. The effects of the (variable) absorption in the 
Earth is taken into account. 

is for this reason that for the detection of astro- 
physical neutrinos it is planned to use an "offline" 
threshold of }Z 1 TeV. The atmospheric back- 
ground above this threshold is small but still po- 
tentially dangerous, it depends on the zenith an- 
gle and is maximum (minimum) for the horizontal 
(vertical) direction at the level of 4 (1) fi/ (km 2 yr). 




Fig. 7. Energy Spectra of ^—induced muons. The thin lines 
are the atmospheric v background integrated in a cone 
of semi-angle 0.3° (dashed lines do not consider oscilla- 
tions). The thick lines are for a point source (a = 2.2, 
$„(> 1 TeV) = lfj- 11 (cm 2 s) _1 . The high (low) line is 
for horizontal (vertical) muons. 

The angular window for the integration of the 
muon signal is determined by three factors: (i) 
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the angular shape of the source, (ii) the intrinsic 
angle #^„, and (iii) the angular resolution of the 
instrument. The source dimension can be impor- 
tant for the galactic sources, in fact the TeV 7-ray 
sources have a finite size, in particular the SNR 
RX J 1713. 7-3946 has a radius - 1°, and Vela Ju- 
nior is twice as large. The detailed morphology of 
these sources measured by the HESS telescope in- 
dicates that most of the emission is coming from 
only some parts of the shell (presumably where the 
gas density is higher), and the detection of v emis- 
sion from these sources could require the careful 
selection of the angular region of the 7-ray signal. 

The distribution of the 8^ v angle defines the min- 
imum angular dimension of a perfect point source 
signal. This is shown in fig. 8 that plots the semi- 
angle of the cone that contains 50, 75 or 90% of 
such a signal (for a a = 2.2 spectrum) as a func- 
tion of the muon energy. The size the muon signal 
shrinks with increasing energy oc (-E^) -1 / 2 . This is 
easily understood noting that the dominant con- 
tribution to d^ v is the muon-neutrino scattering 
angle at the ^-interaction point: 

m N x ( E^ \ 
cos V = 1 - - -g- j (3) 

(E^fi is the muon energy at the interaction point), 
and expanding for small angle one finds 8^ v ~ 
a/tojv x/{2 £7^,0) • The 50% containement cone an- 
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Fig. 8. Angular distribution of the upgoing muon flux in- 
duced by a v power law flux of slope a = 2.2. The curves 
give the semi-angle angle of the cone that (for a fixed value 
of Efj,) contains 50%, 75% and 90% of the signal. The thick 
line gives the cone that maximize the S/^fB. 



gle shrinks to 6» 50% ~ 0.15° for E^ = l TeV that is 
probably smaller than the experimental resolution. 

Qualitatively the experimental strategy for the 
maximization of the point-source sensitivity is 
clear: (i) selection of high energy > 1 TeV 
muons, and (ii) selection of the narrowest angular 
window compatible with the experimental resolu- 
tion and the signal angular size. These cuts will 
reduce the total signal by at least a factor of 2. 
The optimization of the cuts for signal extraction 
is a non-trivial problem, that we cannot discuss 
here. In principle one would like to maximize the 
quantity S/y/B (where S and B are respectively 
the signal and background event numbers). For 
the angular window (where B oc 8) this problem 
has a well defined solution that is shown in fig. 8, 
that corresponds to the choice of a window that 
contains approximately 50% of the signal. This 
is the numerically obtained generalization of the 
well known fact that for a gaussian angular reso- 
lution of width <tq , S/ y[~B is maximized choosing 
the angular window 8 < 1.585 o~e that contains a 
fraction 0.715 of the signal. For the energy cut, 
the quantity S/ y[~B grows monotonically when the 
threshold energy is increased, and therefore the 
optimum is determined by the condition that the 
background approaches zero. 

6.2. Sensitivity and Sources 

The conclusion of this discussion on the point- 
source sensitivity is that after careful work for 
background reduction, the minimum v flux (above 
1 TeV) detectable as a point source for the km 3 
telescopes is of order (or larger than) Q v ~ 
10~ n (cm 2 s) _1 , that is approximately the v flux 
of the reference source considered above. This is 
tantalizingly (and frustratingly) close to the flux 
predicted from the extrapolation (performed un- 
der the assumption of a transparent source) of the 
brightest sources in the TeV 7-ray catalogue. 

The remarkable recent results of HESS and 
MAGIC on galactic sources are certainly very ex- 
citing for the v telescope builders, because they 
show the richness of the "high energy" sky, but 
they are also sobering, because they lead to a 
prediction of the scale of the brightest galactic v 
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sources, that is just at the limit of the sensitivity of 
the planned instruments. As discussed previously, 
strong v sources could lack a bright 7-ray coun- 
terparts because of internal photon absorption, or 
because they are at extra-galactic distances. 

The considerations outlined in this section can 
be immediately applied to a burst-like transient 
sources like GRB's. In the presence of a photon 
"trigger" that allows a time coincidence, the back- 
ground problem disappears and the detectability of 
one burst is only defined by the signal size. A burst 
with an v energy fluence F v ~ 2 x 10 erg/cm 2 
per decade of energy (in the 10-100 TeV range) 
will produce an average of one ^-induced muon in 
a km 3 telescope. 



7. Extragalactic Neutrino Sources 

The flux from the ensemble of all extragalactic 
sources, with the exception of a few sources will ap- 
pear as an unresolved isotropic flux, characterized 
by its energy spectrum <f>(E). The identification of 
this unresolved component can rely on three sig- 
natures: (i) an energy spectrum harder that the 
atmospheric flux, (ii) an isotropic angular distri- 
bution, and (iii) approximately equal fluxes for all 
6 neutrino types. This last point is a consequence 
of space averaged flavor oscillations. 

The flux of prompt (charm decay) atmospheric 
neutrinos is also approximately isotropic in the 
energy range considered, and is characterized by 
equal fluxes of v e and u^, with however a signif- 
icant smaller flux of v T (that are only produced 
in the decay of Df). The disentangling of the as- 
trophysical and an prompt-atmospheric fluxes is 
therefore not trivial, and depends on a good de- 
termination of the neutrino energy spectrum, to- 
gether with a convincing model for charmed par- 
ticle production in hadronic interactions. A model 
independent method for the identification of the 
astrophysical flux requires the separate measure- 
ment of the fluxes of all three neutrino flavors, in- 
cluding the v r . The flux of prompt atmospheric 
is also accompanied by an approximately equal flux 
of n T (the differences at the level of 10% are due to 
the very well understood differences in the spectra 



produced in weak decays). The down-going muon 
flux is in principle measurable, determining the 
prompt component experimentally and eliminat- 
ing a potentially dangerous background. This mea- 
surement is therefore important and efforts should 
be made to perform it. 

7.1. Energetics of Extra-galactic neutrinos 

The most transparent way to discuss the v ex- 
tragalactic flux is to consider its energetics. The 
observable v flux <p(E) is clearly related to a num- 
ber density n(E) — {An /c) <ft(E) and to an energy 
density p{E) = E n(E), that fill uniformly the uni- 
verse. The v energy density at the present epoch 
has been generated by a power source acting dur- 
ing the history of the universe. Integrating over all 
E„, the relation between the power injection den- 
sity C(t) and the v energy density at the present 
epoch is given by: 



to 



dt 



= / dz 



C{t) 



[l + z(t)] 



dz 



C(z) 
(1 + *) 



/ 



C{z) 



Co 



H(z) (1 + zf H 



(4) 



where to is the age of the unverse, Co and Ho are 
the v power density and the Hubble constant at the 
present epoch, and £ is the adimensional quantity: 



dz 



Hp 

H(z) 



Co 



(1 + *)- 



(5) 



that depends on the cosmological parameters, and 
more crucially on the cosmic history of the neu- 
trino injection. The quantity Co £ can be under- 
stood as the effective "average power density" op- 
erating during the Hubble time (.Ho) -1 , and de- 
pends on the cosmological parameters, and more 
strongly on the cosmological history of the injec- 
tion. For a Einstein-De Sitter universe (fi m = 1, 
Qa = 0) with no evolution for the sources one has 
£ = 0.4, for the "concordance model" cosmology 
(fi m = 0.3, tt A = 0.7) this becomes £ = 0.53. 
For the same concordance model cosmology, if it 
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is assumed that the cosmic time dependence of 
the neutrino injection is similar to the one fitted 
to the star-formation history [25,26], one obtains 
£(SFR) ~ 3.0, for a time dependence equal to the 
one fitted to the AGN luminosity evolution [27] one 
finds £(AGN) ~ 2.2. 

If we consider the power and energy density not 
bolomctric, but only integrated above the thresh- 
old energy E m i n , the redshift effects depend on the 
shape of the injection spectrum. For a power law 
of slope a independent from cosmic time, one has 
still a relation of form: 



-no 

very similar to (4) , with: 



oo 

/ 



dz 



\ H - 






[H(z)_ 




I 4) J 



(6) 



(7) 



For a — 2 one has £2 = £• The energy density (6) 
corresponds to the v isotropic differential flux: 



b v {E) = K v E- a 
£0 (E n 



1 ^min £<* k a 



E~ 



(8) 



4ttH 

where k a is the kinematic adimensional factor: 
k a ~ (a - 2)/[l - (E^/E^)*- 2 ]- 1 (9) 



j max the high energy cutoff of the spec- 



with E B _ 

trum. The limit of k a when a — > 2 is: k a = 
[log(-E max /.E m i n )] _ . One can use equation (8) to 
relate the observed diffuse v flux to the (time and 
space) averaged power density of the neutrino 
sources: 



(K, E^) 



4tt H t 







(10) 



The case a = 2 corresponds to an equal power 
emitted per energy decade, and numerically: 



K v ~ 3.7 x 10" 



(£0 C)decade 



L /Mpc J 



GeV 



cm s sr 



(11) 



(Lq is the solar luminosity). The current limit on 
the existence of a diffuse, isotropic v flux of slope 
a ~ 2 from the AMANDA and Baikal detectors 



corresponds to K v < 8 x 10 -7 GeV/(cm 2 s sr) -1 , 
while the sensitivity of the future km 3 telescopes 
is estimated (in the same units) of order K v ~3x 
10~ 9 . This implies that the average power density 
for the creation of high energy neutrinos is: 

(£oO < 2.4 x 10 4 (L e Mpc" 3 ) decade -1 (12) 

the future km 3 telescopes will detect the extra- 
galactic flux if it has been generated with an "av- 
erage" power density (£0 ~ 80 (same units). 



7.2. Possible Sources of Power 

To evaluate the astrophysical significance of 
the current limits on the diffuse u flux and of the 
expected sensitivity of the km 3 telescopes, one 
can consider the energetics of the most plausible 
sources of high energy neutrinos, namely the death 
of massive stars and the activity around Super 
Massive Black Holes at the center of galaxies. 

7.2.1. The Death of Stars 

The bolometric average power density asso- 
ciated with stellar light has been estimated by 
Hauser and Dwek [28] as: 



(£?star) b01 



(0.36 -e- 1.23) x 10° (L Q Mpc -3 ) (13) 



this is also in good agreement with the estimate of 
Fukugita and Peebles [26] for the B-band optical 
luminosity of stars. This power implies that the 
total mass density gone into the star formation [26] 
is: 



Psta 



5.9 x 10 s (Mq Mpc" 



(14) 



(that implies fi sta r = 0.0043 ~ 0.09 fi ba ryon)- Most 
of this power has of course little to do with cosmic 
rays and high energy v production, however c.r. ac- 
celeration has been associated with the final stages 
of the massive stars life. A well established mech- 
anism is acceleration in the spherical blast waves 
of SuperNova explosions, a more speculative one 
is acceleration in relativistic jets that are emitted 
in all or a subclass of the SN, and that is the phe- 
nomenon behind the (long duration) Gamma Ray 
Bursts. The current estimate of the rate of gravita- 
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tional collapse supernova rate obtained averaging 
over different surveys [26] is: 



^observed _ 7 g+6.4 x 1Q -4 ( Mp( ,3 yJ .yl 



This is reasonably consistent with the the estimate 
of the power emitted by stars, if one assumes that 
all stars with M > 8 M end their history with a 
gravitational collapse. Fukugita and Peebles esti- 
mate (V'(O) is the SFR at thre present epoch): 



Ray = 1>(0) 



100 



AM dN_ 

alV1 dM 



100 
0.08 



dMMj§ 



■7.9t 2 3 4 9 x 10" 



(Mpc 3 yr)" 



(16) 



Assuming that the average kinetic energy released 
in a SN explosion is of order (-Ekin)siv — 1-6 x 
10 51 erg this implies an average power density: 



SN,kin 



: 4.2 x 10 6 (L Q Mpc 



(17) 



It is commonly believed that a fraction of order 
~ 10% of this kinetic energy is converted into rel- 
ativistic hadrons. One can see that if on average a 
fraction of order 1% is converted in neutrinos, this 
could result in (£ 0f ~ 10 3 (L Q Mpc -3 ), that 
could give a a detectable signal from the ensemble 
of ordinary galaxies. 



7.2.2. Active Galactic Nuclei 

Active Galactic Nuclei have often been sug- 
gested as a source of Ultra High energy Cosmic 
Rays, and at the same time of high energy neutri- 
nos. Estimates of the present epoch power density 
from AGN in the [2,10 KeV] A-ray band are of 
order 

4 2 ^ cV1 ~ 2.0 x 10 5 (L Mpc -3 ) (18) 

with a cosmic evolution [27] characterized by £ ~ 
2.2. Ueda [27] also estimates the relation between 
this energy band and the bolomctric luminosity is: 
iboi — 30 lJ£' 10 KcV ' . and therefore the total power 
associated with AGN is 

(££)agn^10 7 (L Mpc- 3 ) (19) 

of order ~ 1% of the power associated with star 
light. It is remarkable that this AGN luminosity is 



well matched to the average mass density in Super 
Massive Black Holes, that is estimated [29,30]: 



(15) p. = A.6±\i x 10 5 M Q Mpc" 



(20) 



When a mass m falls into a Black Hole (BH) of 
mass M, the BH mass is increased by an amount: 
AM. = (1 — e) m while the energy e m is radiated 
away in different forms. The total accreted mass 
can therefore be related to the amount of radiated 
energy, if one knows the radiation efficiency e: 



M. 



(1-e) 



E 



radiated 



and therefore 



(£0 



AGN H Q 



(21) 



(22) 



The estimates (19) and (20) are consistent for an 
efficiency s ~ 0.1. Assuming radiated energy is 
the kinetic energy accumulated in the fall down to 
a radius r that is / times the BH Schwarzschild 
radius R s = 2GM.: 



E, 



e = 



radiated 

m 



GM. _ 1 



(23) 



this corresponds to / ~ 5. 

The ensemble of AGN has sufficient power to 
generate a diffuse v flux at the level of the existing 
limit. 

7.3. Gamma Ray Bounds 

A very important guide for the expected flux of 
extragalactic v is the diffuse 7 ray flux measured 
above 100 MeV measured by the EGRET instru- 
ment [8] that can be described (integral flux above 
a minimum energy E 1 in GeV) as: 



$ 7 (£; 7 ) ~ 1.42 x lO-^-^tcn^ssr)- 1 



(24) 



(for a critical view of this interpretation see [32] ) . 
This corresponds to the energy density in the 
decade between E 1 and 10 E 7 (with E 1 in GeV): 
p 7 (£ 7 ) ~ 1.35 x 10~ 6 E^ 1 eV cuT 3 and to the 
average power density per decade: 



oca 



3.7 x 10 4 E. 



-0.1 



7 



(L Q Mpc" 



(25) 



This diffuse 7-ray has been attributed [33] to 
the emission of unresolved Blazars. The sum of 
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the measured fluxes (for E 1 > 100 MeV) for 
all identified EGRET blazars [31] amounts to: 
l^Biazars^j — 0-33 (m 2 s) _1 This has to be com- 
pared with the diffuse flux (24) that can be re- 
expressed as: 4tt <I> 7 (100 MeV) ~ 2.62 (m 2 s)" 1 . 
Other authors [34] arrive to the interesting con- 
clusion that unresolved blazars can account for at 
most 25% of the diffuse flux, and that there must 
be additional sources of GeV photons. 

Most models for the Blazar emission favor lep- 
tonic mechanism, with the photons emitted by the 
Inverse Compton scattering of relativistic electrons 
with radiation fields present in the jet, however it 
is possible that hadrons account for a significant 
fraction (or even most) of the 7-ray flux, and there- 
fore that the v flux is of comparable intensity to 
the 7 flux. The test of the "Egret level" for the dif- 
fuse v flux is one of the most important goals for 
the neutrino telescopes already in operation. 

7.4. Cosmic Rays Bounds 

Waxman and Bahcall [35] have suggested the 
existence of an upper bound for the diffuse flux 
of extra-galactic neutrinos that is valid if the v 
sources are transparent for cosmic rays, based on 
the observed flux of ultra high energy cosmic rays. 
The logic (and the limits) of the WB bound are 
simple to grasp. Neutrinos are produced by cosmic 
ray interactions, and it is very likely (and certainly 
economic) to assume that the c.r. and v sources 
are the same, and that the injection rates for the 
two type of particles are related. The condition of 
"transparency" for the source means that a c.r. 
has a probability P- m t < 0.5 of interacting in its 
way out of the source. The transparency condi- 
tion obviously sets an upper bound on the v flux, 
that can be estimated from a knowledge of the 
cosmic ray extragalactic spectrum, calculating for 
each observed c.r. the spectrum of neutrinos pro- 
duced in the shower generated by the interaction 
of one particle of the same mass and energy, and 
integrating over the c.r. spectrum. If the c.r. spec- 
trum is a power law of form <fi c , T .(E) = K E~ a , 
the "upper bound" v flux obtained saturating the 
transparency condition (for a c.r. flux dominated 
by protons) is: <j) v = (K Z p ^ u ) E~ a with Z v ^, v ~ 



0.25 

The WB upper bound (shown in fig. 9) has been 
the object of several criticisms see for example [36] . 
There are two problems with it. The first one is 
conceptual: the condition of transparency is plausi- 
ble but is not physically necessary. The c.r. sources 
can be very transparent, for example in SNR's the 
interaction probability of the hadrons accelerated 
by the blast wave, is small (<C 1%) (with a value 
that depends on the density of the local ISM), but 
"thick sources" are possible, and have in fact been 
advocated for a long time, the best example is ac- 
celeration in the vicinity of the horizon of a SMBH. 
The search for "thick" neutrino sources is after all 
one of the important motivations for v astronomy. 

The second problem is only quantitative. In or- 
der to estimate the bound one has to know the 
spectrum of extra-galactic cosmic rays. This flux 
is "hidden" behind the foreground of galactic cos- 
mic rays, that have a density enhanced by mag- 
netic confinement effects. The separation of the 




E (eV) 



Fig. 9. Different contributions to the average energy density 
in the universe. Different lines describe the contributions of 
the 2.7° K CMBR radiation, starlight (with reproccssion 
by dust), the X-ray diffuse flux (attributable to AGN, the 
Egret (_E 7 > fOO MeV) diffuse flux, and the GRB contri- 
bution from the BATSE data. The detected cosmic ray en- 
ergy density (points) is enhanced by magnetic confinement 
in our galaxy and is not cosmological, except (possibly) at 
the highest energy. The lines labeled WB(cr) and WB(i/) 
are the extragalactic c.r. and upper bound v fluxes esti- 
mated in [35] . The line labeled Amanda is the current limit 
on an isotropic extragalactic v flux, and the line labeled 
ICECUBE is the predicted sensitivity on km 3 telescopes. 

galactic and extragalactic component of the c.r. is 
a central unsolved problem for cosmic ray science. 
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The extragalactic c.r. is dominant and therefore 
visible only at the very highest energies, perhaps 
only above the "ankle" (at -E an kic ~ 10 19 eV) as 
assumed in [35] . More recently it has been argued 
[37] that extra-galactic p dominate the c.r. flux 
for Eq > 10 18 eV. Waxman and Bahcall have fit- 
ted the c.r. flux above the ankle with an a E~ 2 
injection spectrum (correcting for energy degra- 
dation in the CMBR) and extrapolated the flux 
with the same form to lower energy. The power 
needed to generate this c.r. density was £ cr ~ 
1900 ^©/(Mpc 3 decade). The extrapolation of the 
c.r. flux (that is clearly model dependent) and the 
corresponding v upper bound are shown in fig. 9, 
where they can be compared to the sensitivity of 
the v telescopes. 

Because of the uncertainties in the fitting of the 
c.r. extragalactic component and its extrapolation 
to low energy (and the possible loophole of the ex- 
istence of thick souces) the WB estimate cannot 
really be considered an true upper bound of the 
v flux. However the "WB i/-flux" is important as 
reasonable (and indeed in many senses optimistic) 
estimate of the order of magnitude of the true flux. 
The existence of cosmic rays with energy as large 
as E ~ 10 20 eV is the best motivation for neutrino 
astronomy, since these particles must be acceler- 
ated to relativistic energies somewhere in the uni- 
verse, and unavoidably some c.r. will interact with 
target near (or in) the source producing neutrinos. 

7.5. Resolved and Unresolved Fluxes 

An interesting question is the relation between 
the intensity of the total extragalactic v contribu- 
tion, and the potential to identify extra-galactic 
point sources. This clearly depends on the lumi- 
nosity function and cosmological evolution of the 
sources. Assuming for simplicity that all sources 
have energy spectra of the same shape, and in 
particular power spectra with slope a, then each 
source can be fully described by its distance and 
its v luminosity L (above a fixed energy threshold 
-Emm)- The ensemble of all extra-galactic sources 
is then described by the function n(L, z) dL that 
gives the number of sources with luminosity in the 
interval between L and L + dL contained in a unit 



of comoving volume at the epoch corresponding to 
redshift z. The power density due to the ensemble 
of all sources is given by: 



jC(z) = J dL L n(L, z) 



(26) 



It is is possible (and in fact very likely) that 
most high energy neutrino sources are not be 
isotropic. This case is however contained in our 
discussion if the luminosity L is understood as 
an orientation dependent "isotropic luminosity": 
L = (47r) (dL tIue /dSl). For a random distribution 
of the viewing angles it is easy to show that 



dL L n(L, z) = dL tIUC L truc n(L tIue , z) (27) 



The v flux (above the threshold energy E m [ n ) re- 
ceived from a source described by L and z is: 



point 



(L,z) 



L 



k a {l + zf 



(a-1) d L {zf 



(28) 



where cLl{z) is the luminosity distance and k a is 
the adimensional factor given in (9). 

If ^min is the sensitivity of a neutrino telescope, 
that is the minimum flux for the detection of a 
point source, then a source of luminosity L can be 
detected only if is closer than a maximum distance 
corresponding to redshift Zh(L). Inspecting equa- 
tion (28) it is simple to see that Zh is a function of 
the adimensional ratio x — yjL/L* where: 



L* = 4^£ min $ min (a - 1)/ (Hi k a 



(29) 



is the order of magnitude of the luminosity of a 
source that gives the minimum detectable flux 
when placed at z <~ 1. The explicit solution for 
Zh(x) depends on the cosmological parameters 
(f2 m , Oa) and on the spectral slope a. A general 
closed form analytic solution for z h (x) does not 
exist 6 , but it can be easily obtained numerically. 
The solution for the concordance model cosmol- 
ogy is shown in fig. 10. It is also easy and useful to 
write Zh as a power law expansion in x: 



z h = x+- [2(1 - Q A ) 



0, 



2a]x 2 + ... (30) 



6 As an example, for f2 m = 1, f2j\ = and a = 2 the exact 
solution is = (x — 1 + \/l + 2x)/2. 
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Points: Egret Blazars [0-1-1 GeV] 



Fig. 10. Redshift horizon Zh(L) calculated for the cosmol- 
ogy f2 m = 0.3, Oj\ = 0.7, and sources with power law spec- 
tra of slope a = 2, assuming the a minimum detectable flux 
#i/(> 1 TeV = 10 -11 (cm 2 s) — 1 . The points corresponds 
to the EGRET blazars [31] with the luminosity measured 
in the 0.1—10 GeV range. 



The leading term of this expansion is simply x in- 
dependently from the cosmological parameters and 
the slope a. This clearly reflects the fact that for 
small redshift z one probes only the near universe, 
where and when redshift effects and cosmological 
evolution are negligible, and the flux simply scales 
as the inverse square of the distance. 

The total flux from all (resolved and unresolved) 
sources can be obtained integrating over L and z: 



*tot = ^- / <IL 

47T 



dz v(z) n(L, z) $p i nt (i, z)(31) 



where v(z) is the comoving volume contained be- 
tween redshift z and z + dz: 



v(z) = Air 



dl(z) 



(1 



H(z) 



4-7T 



(! + ■■•) (32) 



Substituting the definitions (28) and (32), integrat- 
ing in L and using (26) the total flux corresponds 
exactly to the result (8) . The resolved (unresolved) 
flux can be obtained changing the limits of inte- 
gration in redshift to the interval z £ [0, Zh(L)] 
(z £ [zh(L),oc]). The total number of detectable 
sources is obtained integrating the source density 
in the comoving volume contained inside the hori- 
zon: 



iVs 



dL 



dz v{z) n(L, z) 



(33) 



A model for n(L, z) to describe the luminosity dis- 
tribution and cosmic evolution of the sources al- 
lows to predict the fraction of the extra-galactic 
associated with the resolved flux, and the corre- 
sponding number of sources. Here there is no space 
for a full discussion, but it may be instructive to 
consider a simple toy model where all sources have 
identical luminosity (that is n{L,z) — n 8[L — 
L(z)]). In this model, is the (unique) luminosity L 
of the sources at the present epoch is not too large, 
(that is for L < L* with L* given in (29)), for a, fixed 
total flux for the ensemble of sources, the number 
of objects that can be resolved is (for a ~ 2): 



iVs 



1 H K v y/L 



decade 



3 £, log(10) c ($ min £ min ) 3 / 2 
* 10 (K_ 9 ) (O'Tir 1 (£«)* y (34) 

where K v is the coefficient of the diffuse v flux 
(K_$ is in units 10 -9 GeV/(cm 2 s sr)), idoc is the 
power of an individual source per energy decade 
(L 45 is in units of 10 45 erg/s), and $ m i n is the 
minum flux above energy E m i n = 1 TeV for source 
identification (^-n is in units 1CP 11 (cm 2 s) _1 ). 
It is easy to understand the scaling laws. When 
the luminosity of the source increases the radius 
of the source horizon (for L not too large) grows 
oc y i/'&min and the corresponding volume grows 
as V cx (i/$min) 3 ^ 2 , while (for a fixed total flux) 
the number density of the sources is n cx K v L , 
therefore the number of detectable sources is: N s cx 
nVcx KvL 1 ' 2 $J n /2 . 

Similarly, the ratio of the resolved to the total v 
flux can be estimated as: 



^resolved 
^total 



1 1 

! 4^ \ log(10) 
• 0.005 ($ 



mm 



Hi 

c 2 
-1 



^dec 



^min 

(T \ 2 - 5 



(35) 
(36) 



The scaling oc L/ <& m i n also easily follows from the 
assumption of an euclidean near universe. 

The bottom line of this discussion, is that it is 
very likely that the ensemble of all extra-galactic 
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sources will give its largest contribution as an un- 
resolved , isotropic contribution, with only a small 
fraction of this total flux resolved in the contribu- 
tion of few individual point sources. The number of 
the detectable extra-galactic point sources will ob- 
viously grow linearly with total extra-galactic flux, 
but also depends crucially on the luminosity func- 
tion and cosmological evolution of the v sources. 
If a reasonable fraction of the individual sources 
are sufficiently powerful (L > 10 45 erg/s), an inter- 
esting number of objects can be detected as point 
sources. Emission for blazars is a speculative but 
very exciting possibility (see fig. 10). 



8. GZK neutrinos 

For lack of space we cannot discuss here the 
neutrinos of "GZK" origin and other speculative 
sources [38]. This field is of great interest, and 
is in many ways complementary to the science 
that can be performed with km 3 telescopes. GZK 
neutrinos (see fig.l) have energy typically in the 
10 18 -10 20 eV, and the predicted fluxes are so small 
that km 3 detectors can only be very marginal, and 
larger detector masses and new detection methods 
are in order. Several interesting ideas are being 
developed (for a review see [39]), these include 
acoustic [40], radio [41] and Air Shower [42,43] 
detection. 

The GZK neutrinos are a guaranteed source, and 
their measurement carry important information on 
the maximum energy and cosmic history of the 
ultra high energy cosmic rays. Perhaps the most 
promising detection technique uses radio detectors 
on balloons. An 18.4 days test flight of the ANITA 
experiment [7] (see fig. 1) has obtained the best 
limits on the diffuse flux of neutrinos at very high 
energy. 

As photon astronomy is articulated in different 
fields according to the range of wavelength ob- 
served, for neutrino astronomy one can already see 
the formation of (at least) two different subficlds: 
the "km 3 neutrino Science" that aims at the study 
of v in the 10 12 -10 16 eV energy range, and "Ul- 
tra High Energy neutrino Science" that studies v 
above 10 18 eV, with the detection of GZK neutri- 



nos as the primary goal. 



9. Conclusions and Outlook 

The potential of the km 3 telescopes to open the 
extraordinarily interesting new window of high en- 
ergy neutrino astronomy is good. The closest thing 
to a guaranteed v source are the young SNR ob- 
served by HESS in TeV photons. The expected 
v fluxes from these sources are probably above 
the sensitivity of a km 3 telescopes in the north- 
ern hemisphere. Other promising galactic v sources 
are the Galactic Center and /iQuasars. Blazars and 
GRB's are also intensely discussed as extragalac- 
tic v sources. The combined v emission of all ex- 
tragalactic sources should also be detectable as a 
diffuse flux distinguishable from the atmospheric 
v foreground. Because of the important astrophys- 
ical uncertainties, the clear observations of astro- 
physical neutrinos in the km 3 telescopes is however 
not fully guaranteed. 

There are currently plans to build two different 
instruments of comparable performances (based on 
the water Cherenkov technique in water and ice) 
in the northern and southern hemisphere. Such in- 
struments allow a complete coverage of the celes- 
tial sphere. This is a very important scientific goal, 
if the detector sensitivity is sufficient to perform 
interesting observations. If the deployment of one 
detector anticipates the second one, it is necessary 
to be prepared to modify the design of the second 
one, on the basis of the lessons received. This is 
particularly important if the first observations give 
no (or marginal) evidence for astrophysical neutri- 
nos, indicating the need of an enlarged acceptance. 

A personal "guess" about the most likely out- 
come for the operation of the km 3 telescopes, is 
that they will play for neutrino astronomy a role 
similar to what the first A-ray rocket of Rossi and 
Giacconi[l,2] played for A-ray astronomy in 1962. 
That first glimpse of the X-ray sky revealed one 
single point source, the AGN Sco-Xl (that a the 
moment was in a high state of activity), and ob- 
tained evidence for an isotropic A-ray light glow 
of the sky. Detectors of higher sensitivity (a fac- 
tor 10 4 improvement in 40 years) soon started to 
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observe a large number of sources belonging to dif- 
ferent classes. 

Even in the most optimistic scenario, the 
planned km 3 telescopes will just "scratch the 
surface" of the rich science that the neutrino mes- 
senger will carry. To explore this field it will be 
obviously necessary to develop higher acceptance 
detectors, and it is not too soon to think in this 
direction. 
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